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Abstract 

According to the general theory of relativity, kinetic energy contributes 
to gravitational mass. Surprisingly, the observational evidence for this 
prediction does not seem to be discussed in the literature. I reanalyze 
existing experimental data to test the equivalence principle for the 
kinetic energy of atomic electrons, and show that fairly strong lim- 
its on possible violations can be obtained. I discuss the relationship 
of this result to the occasional claim that "light falls with twice the 
acceleration of ordinary matter." 
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The principle of equivalence — the exact equality of inertial and gravitational 
mass — is a cornerstone of general relativity, and experimental tests of the univer- 
sality of free fall provide a large set of data that must be explained by any theory 
of gravitation. But the implication that energy contributes to gravitational mass 
can be rather counterintuitive. Students are often willing to accept the idea that 
potential energy has weight — after all, potential energy is a rather mysterious 
quantity to begin with — but many balk at the application to kinetic energy. Can 
it really be true that a hot brick weighs more than a cold brick? 

General relativity offers a definite answer to this question, but the matter is 
ultimately one for experiment. Surprisingly, while observational evidence for the 
equivalence principle has been discussed for a variety of potential energies, the 
literature appears to contain no analysis of kinetic energy. The purpose of this 
paper is to rectify this omission, by reanalyzing existing experimental data to look 
for the "weight" of the kinetic energy of electrons in atoms. I will then try to 
reconcile the results with the occasional (and not completely unreasonable) claim 
that "objects traveling at the speed of light fall with twice the acceleration of 
ordinary matter." 



1 The Equivalence Principle and Internal Energies 

Modern tests of the principle of equivalence begin with variations of Galileo's apoc- 
ryphal free fall experiment, in which the gravitational accelerations of two objects 
with different compositions are compared. The acceleration due to gravity is pro- 
portional to the ratio mc/raj of gravitational to inertial mass, and violations of 
the principle of equivalence would manifest themselves as differences in accelera- 
tion, or "nonuniversality of free fall." Some experiments!!! a u directly measure free 
fall; otherJH use torsion pendulums to compare centrifugal acceleration to the 
acceleration due to gravity. Galileo's pendulum experiments^ may have alreadv 
achieved accuracies of about 2 x 1CP 3 , and since the pioneering work of E6tv6s,H 
modern experiments have pushed uncertainties down to about 10~ 12 . 

These experiments directly test the equality of gravitational and inertial mass 
for a variety of substances. But because different materials have different compo- 
sitions, the experiments also test the principle of equivalence for various forms of 
internal energy. For example, the inertial mass of an iron nucleus is about 1% less 
than the inertial masses of its constituent protons and neutrons, largely because 
of the (negative) contribution of nuclear binding energy. If this binding energy 
did not also affect gravitational mass, the ratio mo /mi for iron would be greater 
than one, and iron would fall more rapidly than, for example, hydrogen. 

A convenient measure of potential violations of the equivalence principle is the 
Eotvos ratio! 

n(A B) ~ mG{A) mG{B) (I) 
V[A ' B) - mi (A) mi (B)' [ ) 

where mc and mi are the gravitational and inertial masses of two materials A 
and B. Typical experimental limits, which we shall use later, are 

?7(Be,Cu) = (-1.9 ±2.5) x 1(T 12 (2) 
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for beryllium and coppeii and 

i](Al, Pt) = (-0.3 ± 0.9) x 10~ 12 (3) 

for aluminum and platinumJl Now, the inertial mass of a material includes a 
number of contributions — rest mass, nuclear binding energy, electrostatic energy, 
the kinetic energy of constituents, etc. If any of these internal energies E a were to 
violate the principle of equivalence, one could write, to lowest order in the energies, 

m G = mi + XN""HT> ( 4 ) 

a c 

where the r] a are dimensionless parameters that measure the strength of the vio- 
lation. From (|l|) and @, we see that 

V (A,B) = y> ( - Al^l) ■ (5) 

In the absence of fortuitous and fine-tuned cancellations among various types of 
energy, observational limits on rj(A, B) can thus be used to obtain limits on the 
parameters r\ a . 

Expression @ is, of course, an oversimplification: in practice, we can rarely 
separate the internal energy of a material so cleanly into its components. The 
binding energy of protons and neutrons in the nucleus, for example, is an important 
part of the energy of an atom, typically amounting to a bit less than 1% of the 
mass. But protons and neutrons are made up of quarks, and according to quantum 
chromodynamics, nuclear binding energy is merely a remnant of the interaction 
energy of quarks and gluons. Should the binding energy in (Q) include only the 
energy of nucleons, or should we add quark interactions? Should we consider the 
kinetic energy of nucleons, or of the quarks and gluons they comprise? What if 
quarks themselves have constituents? 

Clearly, a complete decomposition of the form (Q) would require a complete 
understanding of the physics of the materials we are studying. Nevertheless, we 
can obtain important, if incomplete, information by isolating a few types of energy 
that are relatively well understood, making the assumption that there will be 
no precise cancellations between these contributions and those of the energies 
we ignore. Without understanding the details of quantum chromodynamics, we 
cannot make sweeping statements about the equivalence principle for all binding 
energy; but we can draw conclusions about the particular contribution of the 
binding energy of nucleons in the atomic nucleus, which is understood empirically 
to a rather high accuracy.0 Similarly, when we consider kinetic energy below, we 
shall not attempt to analyze the kinetic energy of quarks, or even nucleons, but 
shall instead focus on the well- understood physics of atomic electrons. 

Note that even for well-understood components of internal energy, the energy 
of a material is often not directly observable, but must be modeled theoretically. 
We cannot, for example, dismantle the Earth to measure its gravitational binding 
energy, or catch a single atomic electron to measure its kinetic energy while it is 
still in the atom. Often, however, there are useful internal checks on computed 
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energies. In particular, the virial theorem provides an important consistency check, 
and can sometimes be used to relate the energies in which we are interested to 
directly observable quantities. Recall that for a nonrelativistic bound system of n 
particles at positions rj, having total kinetic energy T and interacting through a 
potential U(r), the virial theorem states that 

1 n 

(T) = t J2 (r, • VU) , (6) 

z i=i 

where the angle brackets denote a time average. In particular, for an electromag- 
netically bound system, U(r) ~ J2Qi/\ r ~ r «l; an d it is not hard to show that 

T = -U/2 = -E, (7) 

where U is the electrostatic potential energy and E is the total (kinetic plus 
potential) energy. We shall use this relation below to check computed values of 
kinetic energy for atomic electrons. 

Using an analysis based on equation (||) , physicists have obtained strong limits 
on violation of the equivalence principle by strong interaction energy in nuclei, 
electrostatic and magnetostatic energy in nuclei, the energy of hyperfine interac- 
tions, and weak interaction energyEI The^lectrostatic binding energy of electrons 
in atoms has also been briefly discussedJill A somewhat different method, which 
uses Lunar laser ranging to compare the accelerations of the Earth and the Moon, 
has led to limits on violations of the equivalence principle by gravitational binding 
energy.EJ 



2 The Case of Kinetic Energy 

Our interest in this article is kinetic energy, which has, surprisingly, not yet been 
analyzed in the literature. Most of the kinetic energy of an atom resides in the 
nucleus: a typical jmcleus of atomic number A has a radius R ~ RqA 1 ^ 3 with 
Rq ~ 1.3x 10~ 15 m,H3 and the uncertainty principle yields an estimate T nuc /mic 2 ~ 
10 -2 A~ 2 / 3 . This argument only sets a lower bound on the kinetic energy, however, 
and may not give the actual dependence on A; to use equation (|5|), we would 
need a much more sophisticated and model-dependent calculation. For a system 
described by a simple enough potential, the virial theorem relates the kinetic 
energy to the (observable) binding energy. In the nucleus, however, several types of 
potential energy compete, and the virial theorem does not give a unique separation 
of energies. A decomposition like that of equation (|j) is therefore problematic; we 
do not know enough to distinguish the kinetic energy from other contributions. 

For electrons in atoms, on the other hand, these problems largely disappear. 
Accurate computations of electron kinetic energies are now standard in condensed 
matter physics, and while exact solutions of the many-body Schrodinger equation 
are not known, well-understood and well-tested approximations are readily avail- 
able. Moreover, atomic electrons are bound solely by electromagnetic interactions, 
and the virial theorem may be used to check computed kinetic energies against ob- 
served binding energies. This simplicity comes at a price: electron kinetic energy 
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is only a small part of total energy, T e ie C /mi(? ~ 10~ 7 , and experiments are thus 
less sensitive to possible violations of the equivalence principle. Nevertheless, the 
existing data are accurate enough to provide a good test in this relatively clean 
system. 

The basic physics we wish to explore is fairly simple. The electrons in an atom 
with high atomic number are, on average, more tightly bound than those in an 
atom with low atomic number. Their kinetic energy is consequently greater, and 
constitutes a greater proportion of the total energy of the atom. If, as an extreme 
example, this kinetic energy had no weight, a high-Z atom would fall measurably 
more slowly than a low-Z atom. 

To obtain quantitative predictions, we need to determine the kinetic energy of 
atomic electrons^ We can begin with the simple Thomas-Fermi model for many- 
electron atoms which treats the electrons statistically as a Fermi gas and uses 
semiclassical methods to determine their characteristics. In this approximation, 
a typical electron in an atom of atomic number Z is located roug hly Z^ 1 / 3 Bohr 
radii from the nucleus. The electrostatic energy of one such electron is proportional 
to Z/Z -1 / 3 = Z 4 / 3 ; the energy of Z electrons thus goes as Z 7 / 3 . By the virial 
theorem (Q), the kinetic energy should have the same form, T ~ cZ 7//3 . From the 
observed ground state energy of the hydrogen atom, we can estimate c to be on 
the order of 10 eV. The ratios T e i ec /mjc 2 in equation (||) thus range from about 
10 -6 for platinum to 10~ 8 for beryllium, much larger than the limits (@)-(0) for 
violation of the equivalence principle. 

To find more precise results, we could numerically integrate the Thomas- Fermi 
model to determine the coefficient c. But it is almost as easy to employ the 
much more accurate numerical approaches that are now widely used in condensed 
matter physics. Consider, for example, the kinetic energies of atomic electrons in 
beryllium and copper. These can be computed in the local density approximation^ 
(for a review of this method, see reference |15|), using standard and widely available 
computer codes. For isolated atoms, one finds that 

-^l(Be) = 4.6x 1(T 8 



mjc 
mic 2 



(Cu) = 7.7 x 1(T 7 . (8) 



To check these numbers, we can appeal to the virial theorem (0) for electro- 
magnetically bound systems, which allows us to compare the energies in equa- 
tion (|8|) to published values of total energies as computed in the Hartree-Fock 
approximation!!!! and the local density approximation. O The results agree to within 
2%. Better yet, the total energy E can be measured directly — it is the ionization 
potential, the energy required to totally ionize an atom — and the kinetic energies 
(|D can be compared to these observations. We again obtain agreement to within 
3% for berylliumta and 1% for copper .Ej One might worry that our computations 
were performed for isolated atoms, while the experimental test of reference ^ used 
solid metallic beryllium and copper. But the relevant energy differences, the co- 
hesive energies of the metals, are only a few eV per atom, completely negligible 
for our purposes. 
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We can now combine equations (|2|), (||), and (|8|), assuming as usual that no 
precise cancellations occur among different forms of energy. We obtain a limit 

\rj T \ < 6 x 1(T 6 (9) 

for violation of the equivalence principle by the kinetic energy of electrons. 

A stronger, although less theoretically certain, limit can be obtained by com- 
paring aluminum and platinum. The local density approximation now gives 

^l(Al) = 2.6 x 1(T 7 
777/ cr 

-^(Pt) = 3.3 x 1(T 6 . (10) 

777/ tr 

For aluminum, these figures are again in good agreement with the theoretical!^' 
and experimentalESI results for total energy. For platinum, however, no experi- 
mentalresults appear to be available, and the Hartree-Fock expression for total 
energytB differs from the local density approximation for kinetic energy by about 
15%. The difference may indicate a problem with the computed value; numer- 
ical errors in the code used for these calculations are likely to be significant for 
core electrons in high-Z atoms. In the absence of a better estimate of theoretical 
uncertainties, let us double this difference, and assume conservatively that the 
calculation ( |i"o| ) for platinum is accurate to within 30%. We can then combine 
equations @, (||), and (|l0|) to obtain 

\rfr\ < 6 x 10~ 7 . (11) 

The limits (g) and (|l|) are several orders of magnitude weaker than the corre- 
sponding results for nuclear energies. Nevertheless, they are surprisingly strong, 
and may be counted as good evidence that the equivalence principle holds for 
kinetic energy. 



3 General Covariance and the Weight of Light 

The results of the preceding section will come as no surprise to experts in relativity. 
But perhaps they should. We have another way of "weighing" kinetic energy: we 
can send a beam of particles past a large mass (the Sun, say) and see how it is 
deflected. It is well known that the deflection of light is twice that predicted by 
Newtonian theory; in this sense, at least, light falls with twice the acceleration of 
ordinary "slow" matter. 

Indeed, the general relativistic deflection for a test narticle with an arbitrary 
velocity v and a large enough impact parameter b is0'§l 



2GM 




(12) 

The corresponding angle in Newtonian gravity depends on the ratio of gravita- 
tional to inertial mass, and it is easy to check that equation ( |12D is just the 
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Newtonian result for a particle with an inertial mass mj and a gravitational mass 



I v 2 \ T 
ma = mi I 1 + I = mj + 2~2. (13) 

For light, the kinetic energy T in this expression should be replaced by the elec- 
tromagnetic energy U, which can be loosely interpreted as the "kinetic energy" of 
photons. 

This argument does not seem to be widely published; instead, many texts rely 
on a simple "Einstein elevator" analysis that actually gives only half the correct 
deflection.Ei But it is not uncommon among students, and appears frequently in 
Internet discussions. In such a simple form, the argument is easily addressed: it 
is attempting to impose Newtonian categories on general relativity, ignoring in 
particular the curvature of space.E-3 To obtain a value for the deflection of light by 
the Sun we must at least implicitly compare measurements of direction in widely 
separated regions of space, and to perform such a comparison correctly, we need 
to take into account the curvature of space between these two regions£3 

There is a slightlv more sophisticated version of this argument that is harder 
to dismiss, however.Ej Rather than sending a beam of light past the Sun, let us 
imagine confining the beam to a mirrored box near the Sun, thus avoiding the 
problem of comparing distant frames of reference. If the deflection (|l2|) truly re- 
flects the "weight of kinetic energy," a light beam with energy U should contribute 
an amount 2U to the gravitational mass of the box. 

We can analyze this situation in the weak field approximation to general 
relativityJHH In this limit, the metric is g^ u ~ t]^ u + ha V , with 

^oo = 20, hij = 2<f)5ij, (14) 

where <f> is the Newtonian gravitational potential. The gravitational coupling to a 
test body with a stress-energy tensor T^ v is thusjj] 

i J h^T^d 3 x « J <j> (r 00 + <%T^') d 3 x. (15) 

For a slowly moving particle with rest mass m and kinetic energy T, we have 

T 00 d 3 x w mc 2 + T, J 5 lj T ij d 3 x w 2T (16) 

to lowest order in velocity.il For an electromagnetic field with energy U, on the 
other hand — or for a "photon" of energy U moving at light speed — the stress- 
energy tensor is traceless, so§l 

-i00j3^. ^ tt I X..rni3jZ 



T w d 6 x kU, J SijTH^x w U. (17) 

These expressions may be checked by considering a gas of particles (or photons) 
in a volume V, for which the diagonal spatial components T n of the stress-energy 

^The coupling (|15|) can be obtained from the action S[i/j,g] for matter in curved spacetime by 



expanding the metric around its flat Minkowski value and noting that SS = gT^Sg^u. 



7 



tensor are equal to the pressure p. In that case, (|T^)-(^) may be recognized as 
the standard result that 

pV = (7 - 1)E (18) 

where E is the total (kinetic) energy and the coefficient 7 is 5/3 for nonrelativistic 
particles, 4/3 for the extreme relativistic case. 

Our "box of light" consists of slowly moving walls with energies of the form 
(|16|) and a beam of light with energy of the form (17). For this system — or for 



any other system consisting of electromagnetic fields and matter — we thus have 

J T 00 d 3 x ^mc 2 + T + U, J 5 lj T ij d 3 x &2T + U. (19) 

If <f> is nearly constant over the system, equation (|l^) can thus be interpreted as 
a coupling of the Newtonian potential to the combination mc 2 + 3T + 2U . For 
a free beam of light, the first two terms are absent, and this analysis yields a 
gravitational deflection of twice the Newtonian value, as desired. But comparing 
section |], we seem to have found Eotvos parameters rjr = 2 and rju = 1, in gross 
violation of the equivalence principle. 

Of course, the coupling (|l9|) does not really lead to a disagreement with ex- 
periment: we are saved by the virial theorem. Equation (0) was derived for a 
system of nonrelativistic charged particles coupled by Coulomb interactions. But 
the relativistic derivation of reference ^9] shows that the same relation holds for 
an arbitrary spatially bounded system of electromagnetically interacting particles 
and electromagnetic radiation, as long as U is now understood to be the total 
electromagnetic energy. For our "box of light," 2T + U therefore vanishes, and 
3T + 2U = T + U = E. The apparent violation of the equivalence principle has 
thus rather mysteriously disappeared. 

Such an exact cancellation should have a fundamental explanation. In general 
relativity, it is a consequence of general covariance. Consider a small electromag- 
netically bound "test body" in a gravitational field, with center of mass x % . Let 
us start with the expression (|l5[) and perform a coordinate transformation 



x i -> s*' = x i + <f>£(x) {x L -x { ), (20) 

where £ is any function that is constant inside the test body and falls rapidly to 
zero outside. In the weak field approximation, the fields h^ v transform asEZI 

hfa, -> hfa, - •q lxp d v 5x p - r) up d^5x p , (21) 

and it is easy to check that the coupling ( p!q ) becomes 

4>[mc 2 + T' + U' + (1 - £(0))(2T' + U')\ + higher order terms, (22) 

where T' and U' are the kinetic and electromagnetic energies in the new coordinate 
system. The coupling to 2T+U can thus be altered arbitrarily, and indeed "gauged 
away," by a change of coordinates inside the test body. This argument is closely 
related to our earlier appeal to the virial theorem: the relativistic virial theorem 
can be derived from conservation of the stress-energy tensor fB which is in turn a 
consequence of general covariance in the weak field limit. 
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Put another way, we have learned that the determination of gravitational mass 
from constituent energies is not entirely coordinate- independent. But the physical 
interaction with gravity cannot depend on a choice of coordinates: the coordinate- 
dependent part of the interaction must vanish. That it does is guaranteed by the 
virial theorem, and indeed, this argument can be viewed as a derivation of the 
virial theorem. We have also seen that within the framework of general relativity, 
coordinates can be chosen to show explicitly that the inertial and gravitational 
masses of a bound system are equal. 

We focused so far on electromagnetic energy. But the coupling (|2^) illustrates 
a general ambiguity in determining the parameters r/ a . In any bound system, 
the virial theorem will require that some linear combination of energies vanish, 
and that some linear combination of the rj a therefore be unmeasurable. For an 
electromagnetically bound system, for instance, a term in equation (|j) of the form 
£(2T + U)/c 2 is inherently unobservable. 

How can we deal with this ambiguity? One answer is to compare systems with 
different interactions and different internal energies. For any single system, the 
virial theorem will give a relationship among energies, and one combination c cJ]a 
will be undetermined. But the appropriate combination will differ from system 
to system. Moreover, for systems with interactions described by nonpolynomial 
potentials, it is evident from equation @ that the coefficients c a will depend 
not only on the internal energies, but on the dynamics — the average values of 
the positions — as well. Observations of nuclei, for example, cannot detect a 
gravitational coupling to the combination 2T + U e i ec + cU nuc [ ear , where c depends 
on mean distances between nucleons. 

If the coefficients rj a have any universal significance, we can now combine the 
limits coming from the energy content of nuclei with those coming from atomic 
electrons to obtain information about t]t alone. In particular, the electrostatic 
Eotvos coefficient rju in nuclei satisfies! 

\ VU \ < 4 x 1(T 10 , (23) 

up to an ambiguous coupling that is different from the ambiguous coupling for 
electrons. It seems safe to assume that there should be no perverse cancellation 
between the gravitational couplings of, say, electron kinetic energy in beryllium 
and nuclear binding energy in platinum. Equation ( |23"|) can thus be combined with 
our earlier analysis to remove any "virial theorem ambiguities" in the limits @ 
and (|TTJ) for kinetic energy. 

We can thus tell our students with confidence that kinetic energy has weight, 
not just as a theoretical expectation, but as an experimental fact. 
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